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Abstract 

Remaining within the pure hydrodynamic approach, we formulate a self-consistent model for simulating the dynamic 
behavior of matter passing through metastable states in the two-phase liquid-vapor region of the phase diagram. 
The model is based on the local criterion of explosive boiling, derived by applying the theory of homogeneous bub- 
ble nucleation in superheated liquids. Practical application of the proposed model is illustrated with hydrodynamic 
simulations of a volumetrically uniformly heated planar layer of fused silica SiOi- Implications for experimentally 
measurable quantities are briefly discussed. A newly developed equation of state, based on the well known QEOS 
model and capable of handling homogeneous mixtures of elements, was used in the numerical simulations. 
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1. Introduction 

Most experiments with warm dense matter are inti- 
mately coupled to an adequate analysis of the hydro- 
dynamic motion of the studied material, and the so- 
lution of the hydrodynamic equations requires knowl- 
edge of the thermodynamic equation of state (EOS). 
In many cases — unless the external energy deposition 
goes rapidly well beyond the latent heat of vaporiza- 
tion — the heated (or cooled) material passes through 
the two-phase region of the phase diagram, where the 
liquid and the vapor phases can coexist (for simplic- 
ity, we do not consider the solid phase here). Part of 
this region (between the binodal and the spinodal) cor- 
responds to metastable states, where the thermodynam- 
ically consistent EOS, based on the reversibility consid- 
erations, is essentially double-valued and poses the fol- 
lowing dilemma for any hydrodynamic simulation: one 
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can either follow the metastable branch (the MS-EOS) 
of the superheated liquid (supercooled vapor), or use the 
fully equilibrium EOS (the EQ-EOS) obtained by the 
Maxwell construction across the van der Waals loops. 
In this paper, by using a simple example of a uniformly 
heated planar foil, we show how this dilemma can be re- 
solved within the framework of hydrodynamics, without 
invoking much more complex and less clearly defined 
equations for the kinetics of phase transition. 

Our basic proposition is to follow the evolution of the 
heated liquid into the two-phase region by applying the 
MS-EOS until the criterion of explosive boiling is ful- 
filled. Once the boiling point is reached, we assume that 
the considered fluid element undergoes a practically in- 
stantaneous and irreversible transition to the EQ-EOS. 
The criterion of explosive boiling (|3H . which plays the 
central role in our model, is derived by applying the the- 
ory of homogeneous bubble nucleation in superheated 
liquids Ht. 

The proposed method is applied to a thin planar layer 
of fused silica (SiOa) that is volumetrically uniformly 
heated at a rate of 10 n J g -1 s _1 on a time scale of 
100 ns. Earlier this configuration was proposed |]2| to 
be used in warm-dense-matter experiments under irradi- 
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ation by intense beams of heavy ions at GSfl Initially, 
in its solid-liquid state, the SiC>2 foil is so thin that its 
thermal expansion proceeds in a quasi-static manner, so 
that the pressure across the foil rapidly relaxes to its 
boundary value. For illustrative purposes, we assumed 
in simulations a boundary pressure of po = 500 bar. But 
even if in real experiments much lower boundary pres- 
sures can be anticipated, our basic conclusions should 
remain valid because this value of po is already quite 
small compared to the critical pressure. 

Although in this study only a very simple and per- 
fectly uniform target configuration has been analyzed, 
the proposed method of passing through the two-phase 
region can be applied to a much broader spectrum of 
experimental conditions. How well suited it turns out to 
be for modeling these conditions, remains to be clarified 
by future work. 

This paper is organized as follows. In section [2] the 
quasi-static regime of expansion of a volumetric ally 
heated planar foil is analyzed in detail. In section [3] 
possible role of quasi-elastic oscillations is discussed. 
The criterion for explosive boiling is derived and dis- 
cussed in detail in section [4] The results of hydrody- 
namic simulations, which illustrate practical application 
of the proposed method for treatment of the two-phase 
region, are presented in section All the hydrodynamic 
simulations have been performed by using a newly de- 
veloped EOS package FEOS, whose brief description is 



given in Appendix A 



2. Quasi-static thermal expansion of a planar foil 

Suppose that an initially uniform planar foil at a den- 
sity po is subject to an external volumetric heating with 
a prescribed energy deposition rate q (energy per unit 
mass per unit time). We assume for simplicity that q 
may depend on time t (starting from t — 0) but not on 
the spatial coordinate x, measured perpendicular to the 
foil surface. The ambient (boundary) pressure is fixed 
at some positive value po > 0, which is assumed to be 
small compared to the critical pressure of the liquid- 
vapor phase transition. In the course of heating the foil 
undergoes expansion, and its thickness / = 1( f) increases 
in time starting from the initial value Iq = 1(0). 

Clearly, if the foil is heated slowly by comparison 
with the sonic time 
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where c s is the speed of sound, it expands quasi- 
isobarically in a quasi-static manner. The solution of 
the hydrodynamic equations, describing the quasi-static 
regime of expansion, is readily obtained by successive 
approximations. 

In the zeroth approximation pressure p is constant in 
space and time 



pit, x) = p(p, T) = po, 



(2) 



while the heating rate q determines temporal variation 
of the specific enthalpy h according to the basic thermo- 
dynamic law 

dh 

■ q, (3) 



dt 



which yields 



hit, x) = hit) = hip, T) = h + J qir) dr. (4) 



From the two equations (0 and (|4]i one knows the zero- 
order values of all thermodynamic variables at any given 
time t > 0; their spatial distribution remains perfectly 
uniform. The fluid velocity u in the zeroth approxima- 
tion is zero. 

Substituting the zero-order fluid density p = pit) into 
the mass conservation law 



dp djpu) = 



(5) 



dt dx 

we obtain the first-order solution for the fluid velocity 
u(t, x) = u q (t, x) = ax, (6) 

where 

d In p d\s\v d\s\l 
a = ait) = = — — = — — ; (7) 



dt 



dt dt 



here and below v = p~ l is the specific volume. In Eq. (0 
we assumed that the value x = corresponds to the 
motionless center of mass, i.e. to the middle of the foil. 

By virtue of Eq. ©, the expansion rate a can be cast 
in a form 



a — —q 



dh 
d\np 



ci 



(8) 



where the dimensionless coefficient T — which is just a 
thermophysical parameter of the studied material — is 
defined as 



iw(-*-V 



(9) 
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One readily verifies that for a perfect gas F - y — 1, 
where y is the adiabat index. For a two-term EOS of the 
form 

e(p, T) = e c (p) + c v T, 

! T\ 2 d€ c ^ r T (10) 

pip, T)=p — + T c pc v T, 
dp 

often applied to describe solids and liquids at high en- 
ergy densities J^t], our F coefficient coincides with the 
Griineisen parameter Yq — provided that the latter is 
assumed to be constant; in Eq. ( TTOt e = h - p/p is the 
mass-specific internal energy, cy is a constant heat ca- 
pacity at fixed volume. 

Fhe physical meaning of the expansion rate a is 
straightforward: a is the time scale over which a 
quasi-isobaric heating results in a significant change of 
the density p. Then, the obvious necessary condition for 
applicability of the quasi-static solution @, (HJi and (0 
is 

at s <k 1. (11) 

In other words, the dimensionless product at s is the 
small parameter which defines the hierarchy of suc- 
cessive approximations in our solution of the hydrody- 
namic equations. From Eq. © we find the quasi-static 
expansion velocity u q i of the foil boundary at x = 1/2 



1 c s 
u q i = -al = — iat s ), 



(12) 



which, as might be expected, is a first-order small quan- 
tity compared to the sound speed c s . 

Inequality dTTb is certainly the necessary but not al- 
ways a sufficient condition for practical application of 
the quasi-static solution. Having substituted the first- 
order velocity © into the hydrodynamic equation of 
motion 

du du _j dp 
dx 



, + li- 
St dx 



P 



0. 



(13) 



one finds that this solution is generally characterized by 
small spatial pressure non-uniformities of the order 



Sp 



■P u % 



pc 2 s iat s ) 2 . 



(14) 



Because the mass-specific heating rate q is assumed 
to be perfectly uniform, the same refers to the mass- 
specific entropy s at all times, and the resulting spatial 
non-uniformity of the density can be evaluated as 



Now, if a particular experiment aims at measuring the 
density difference 



Ap = pit + At) - pit) =a p a At 



(16) 



between two quasi-uniform states, separated by a time 
interval At, the measured difference must be large in 
comparison with the inherent spatial non-uniformity, 
i.e. Ap » dp, which implies 



At » t s iat s ). 



(17) 



From practical point of view condition (TTTT i is not par- 
ticularly restrictive because it allows application of the 
quasi-static solution even on time scales At < t s shorter 
than the sonic time — provided, of course, that no sig- 
nificant quasi-elastic oscillations have been excited (see 
next section). 

In practice, the basic condition ( fTTT i for quasi-static 
approximation can be satisfied with relatively high heat- 
ing rates q ^ 10 11 J g~' s _1 , typical for irradiation by 
heavy ion beams at GSI, if a sufficiently thin foil is used. 
In a typical solid the above value of q implies an expan- 
sion rate of a =* 10 7 s _1 , and the product at s falls in the 
range at s < 0.03 for foils with / < 10 pm. 

If a quasi-statically heated foil remains in a liquid 
state and preserves a sharp surface, its expansion veloc- 
ity u q i (which even under quasi-static conditions may 
amount to 10-100 m s~') is an easily measurable quan- 
tity, which provides direct information about the den- 
sity p at any given time t » t s . Once the energy de- 
position fqdt is known and the surface temperature 
T has been measured, one obtains the thermal, a p = 
-p~' idp/8T) p , and the caloric, a* = -p~' idp/dh) p 
expansion coefficients together with the heat capacity 
c p = idh/dT) p at constant pressure. 

A more advanced experimental setup, consisting of 
several foils, is described in Ref. fl. The general goal is 
to employ volumetrically heated samples for EOS mea- 
surements near the boiling curve fl. In particular, the 
subject of non-congruent phase transitions can be ad- 
dressed. Non-congruence means coexistence of two (or 
more) phases with different chemical composition (sto- 
ichiometry). One of the reasons why S1O2 has been 
proposed as a first candidate for such experiments is 
that it is expected to manifest non-congruent evapora- 
tion Jl 



Thus, the quasi-static solution ©, © maintains 
spatial uniformity down to the second order in at s . 



3. Quasi-elastic oscillations and initialization of the 
hydro simulation 

An obvious problem with the quasi-static solution 
©, (|4j, (O is that one cannot in general satisfy the con- 
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Figure 1: Temporal variation of the surface velocity u; of a 10 fim 
thick Si02 planar foil by a sudden turnon of the volumetric heating 
at a rate q = q m = 10" J g - ' s _1 . The amplitude of quasi-elastic 
oscillation is u osc = 3528 cm/s. 




v (cm7g) 

Figure 2: Equation of state of SiC>2 according to the FEOS model on 
the volume-pressure phase plane. Shown are the isotherms with (the 
EQ-EOS) and without (the MS-EOS) the Maxwell construction; the 
binodal and the spinodal touch each other at the critical point CP. 



dition of zero initial velocity u(0, x) = 0. If, for exam- 
ple, the maximum heating rate is turned on at once at 
t = 0, 

0, t < 0, 
q m , t > 0, 



9(0 = 



(18) 



Eq. (O does not allow to start from zero initial velocity 
— even when q m is relatively small and condition (fTTT > 
is satisfied. The reason lies in the fact that Eq. © is not 
the general but a particular first-order solution of the 
hydrodynamic equations for small-amplitude motions 
with m/c s ~ at s « 1. To obtain a general solution, one 
has to take into account quasi-elastic oscillations, repre- 
sented by an infinite series of eigenmodes for standing 
acoustic waves; see . 



Appendix B 



Fig.[T]shows the time variation of the surface velocity 
U[ by a sudden turnon of the volumetric heating as cal- 
culated numerically with the CAVEAT code. Initially, 
the 10 pm thick SiC>2 foil was at rest. One clearly sees 
that periodic oscillations with an amplitude u osc are su- 
perimposed on a slow "secular" expansion of the foil. 



It can be rigorously proven (see Appendix B i that at 



pressure becomes negative in the contraction phase. 

In this work we deal with a situation where quasi- 
elastic oscillations are unwanted because they only in- 
troduce an unnecessary complication into our illustra- 
tive simulations discussed below. Experimentally, os- 
cillations can be suppressed by raising the heating rate 
q smoothly from zero to its maximum level over a time 
scale of At a~ l . In theoretical analysis they can be 
eliminated by simply assigning the linear initial veloc- 
ity profile (0, predicted by the quasi-static solution, as 
it has been done in the present work. 



4. Criterion for explosive boiling 

4.1. General considerations 

Here we explain our principal idea on how to treat 
the hydrodynamics of metastable states on the exam- 
ple of a thin planar foil of S1O2. The realistic EOS of 
S1O2, used in the hydrodynamic simulations, is shown 
in Fig. [2] in standard representation on the (v, p) plane 



t <k a such oscillations have a period t„ = 21q/c s and and described in detail in Appendix A It has all the 
an amplitude u osc = u q i = \alo. 

Clearly, for some experimental tasks quasi-elastic os- 
cillations might be useful (to infer, for example, the 
sound speed c s from the measured period t of oscil- 
lations), for others — highly undesirable because they 
introduce much larger (first-order in at s ) spatial varia- 
tions of pressure and density than are inherent in the 
quasi-static solution. In particular, for sufficiently low 
boundary pressures po < p$c 2 s {at s ), the interior fluid 



characteristic features of the van der Waals EOS: the 
boiling and the condensation curves make up the bin- 
odal; the region of absolute thermodynamic instability 
is delimited by the spinodal, which lies below the bin- 
odal and touches it at the critical point. 

In our example the initially uniform SiOi foil enters 
the metastable region from the liquid side in a quasi- 
static manner due to a uniform volumetric heating. Ap- 
plication of the quasi-static approximation simplifies the 
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Figure 3: Schematic representation of the evolution of the thermo- 
dynamic state of heated matter on the (v,p) phase plane. B\ — > Br. 
transition to the phase equilibrium. 



analysis because, at least initially, the entire foil is rep- 
resented by a single point on the thermodynamic dia- 
grams. It should be noted however that the proposed 
algorithm for traversing the double-valued EOS region 
has general applicability for any dynamic non-uniform 
fluid, entering the metastable region from either the liq- 
uid or the vapor sides. 

Fig. [3] gives a schematic view of the trajectory of the 
heated foil material on the (v, p) plane. We can assume 
that it starts from point O on the cold curve T — from 
the minimum possible specific volume at a given pres- 
sure p = po. Having crossed the liquid binodal at point 
Bo, it continues into the metastable region along the liq- 
uid MS-EOS branch because at point Bo the boiling time 
is formally still infinitely large. 

As the (v, p) state of the heated foil penetrates deeper 
and deeper into the metastable region of the superheated 
liquid, the thermodynamic barrier for transition into a 
fully equilibrium state, prescribed by the EQ-EOS, is 
lowered and, as a consequence, the time scale for tran- 
sition into a new state very rapidly decreases by many 
orders of magnitude. As a result, at a certain point B\ 
one observes very rapid (explosive) boiling of a super- 
heated liquid from inside, irrespective of the boundary 
conditions, i.e. even in an infinitely large volume of a 
perfectly uniform and pure liquid [JJ . 

Keeping this qualitative picture in mind, our basic 
proposition is as follows: the thermodynamic evolution 



of every fluid element is followed along the MS-EOS 
until the criterion for explosive boiling (rapid conden- 
sation) is fulfilled at a certain point B\ between the bin- 
odal and the spinodal; at point B\ one makes an instan- 
taneous transition from the MS-EOS to the EQ-EOS (to 
the point Bo_ in Fig.O at constant values of the fluid den- 
sity p — v and the mass-specific internal energy e. It 
is presumed that the criterion itself is local, i.e. does not 
depend on the boundary conditions and the thermody- 
namic state of adjacent fluid elements. 

Our rationale for keeping precisely the density and 
the specific internal energy constant during the MS — > 
EQ transition is based on the fact that the time scales for 
the p and e change in the considered class of phenom- 
ena are governed by the hydrodynamics and the rate of 
external heating (cooling). Both these time scales are 
usually much longer than that of the explosive boiling 
(if not at the transition point itself, then immediately 
beyond it at a slightly higher temperature; see Fig. |4] 
below). For comparison, near the spinodal the rate of 
transition to the equilibrium state can roughly be esti- 
mated as a few times the mean intermolecular distance 
divided by the mean thermal velocity of molecules, i.e. 
around 10 12 s — which is practically instantaneous by 
any hydrodynamic measure. 

Fixed values of p and e imply that the B\ — > B2 tran- 
sition will be accompanied by an instantaneous jump in 
the values of pressure, temperature and specific entropy. 
According to the general law of increase of entropy Q, 
the final state Bo must always have a higher specific en- 
tropy than the initial state B\ because any fluid element 
at fixed p and e can be considered as a closed thermody- 
namic system undergoing a transition from a metastable 
state to a more probable fully equilibrium state. Less 
obvious is the direction of the temperature change. We 
expect that normally the temperature should also jump 
to a higher value To > T\ (as it does in the FEOS model 
used in this work), although we have managed to prove 
this fact rigorously only for the analytical van der Waals 



EOS itself (see Appendix C 1 



4.2. Mathematical formulation 

To express the criterion for explosive boiling in math- 
ematical terms, we apply the theory of homogeneous 
bubble nucleation in superheated liquids [£T], [sj] . Accord- 
ing to this theory, thermodynamic fluctuations in a su- 
perheated liquid lead to a spontaneous creation of vapor 
bubbles with different radii r. Analysis of the mechan- 
ical and chemical equilibrium of these bubbles reveals 
that there is a critical value r c of the bubble radius, de- 
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fined by the condition 



but quite adequate for our case, will be [HI 



2a- 
P+-, 



(19) 



where cr is the surface tension, and p e = p e (T) is the sat- 
urated vapor pressure at temperature T (see Fig.0; here 
and below the thermodynamic parameters p, T, p and e 
without indices refer to the superheated liquid obeying 
the MS-EOS. Small bubbles with r < r c are dominated 
by surface tension and rapidly collapse, giving no con- 
tribution to the process of boiling. Bubbles with r > r c 
grow and add to the volume fraction occupied by the va- 
por. Hence, the boiling dynamics is governed by the rate 
J [cm -3 s -1 ] of spontaneous creation of vapor bubbles 
with the critical radius r c . 

By its physical meaning, the nucleation rate J defines 
the time scale of boiling of a superheated liquid: if we 
fix a certain fluid volume V and wait for time t, then 
on average JVt critical vapor bubbles will be sponta- 
neously created in this volume. Hence, if one knows 
that in full equilibrium the vapor phase occupies a vol- 
ume fraction f v and that, in addition, the subsequent 
growth of newly created critical bubbles can be ignored, 
one readily calculates that under stationary conditions 
the transition from the metastable to the absolutely sta- 
ble state occurs after a waiting time 



(20) 



where V c = 4nr\j3 is the volume of a critical bubble. 
Below we show that in the dynamic situations analyzed 
in the present work, where the evolution of the ther- 
modynamic state of every fluid element is governed by 
a combination of hydrodynamics and external heating 
(cooling), the growth of supercritical bubbles can in- 
deed be neglected, and that the effective waiting time 
is given by 



dt\ T ) 



(21) 



The key quantity for evaluating the nucleation rate J 
is the minimum work 



m 47T 2 167TCT 3 

3 3(p e - PK 



(22) 



needed to create a critical bubble. Eqs. ([19} and (l22t 
ignore the small Poynting correction to p e -p, discussed 
in Ref. [8]. By combining the arguments of physical 
kinetics and thermodynamics, different versions of the 
basic formula for J have been derived and investigated 
by different authors []]], of which probably the simplest, 



J=N U M~h 



(23) 



where N = p/m is the number of liquid molecules per 
unit volume, and m is the molecular mass. 

We determine the moment t of the MS — > EQ transi- 
tion from the condition that, when a fluid element jumps 
from a metastable liquid state (p, T) - (p,T\) to an 
equilibrium two-phase state (p, T 2 ) (see Fig.|3}, the total 
fractional volume of all vapor bubbles with supercritical 
radii, accumulated in the fluid by that moment, reaches 
its equilibrium value 



< & = &(p, T) 



Pi(T 2 ) -p 
Pl(T 2 ) - Pv(T 2 ) 



< 1, 



(24) 



prescribed by the EQ-EOS for the final state (p, T2); 
here pi(T) and p v (T) are, respectively, the liquid and 
the vapor densities along the corresponding branches of 
the binodal (see Fig. [3). The post-transition temperature 
T2 = T 2 {p, T) is implicitly defined by the equation 



e(p, T) = € E q(p, T 2 ). 



(25) 



In other words, for a given EOS of the van der Waals 
type, T 2 (p, T) is a known function of the density p and 
temperature T in the metastable state. Normally one 
has T 2 (p, T) > T\ for the van der Waals EOS itself the 
mathematical proof of this fact is given in |Appendix C| 
After we take into account that a bubble, created at a 
time f < t with a critical volume V c {t') = 4nrl(f)/3, by 
the time t will occupy the volume 



V(t) 



= V c (t') + j V c (t',t")dt", 



(26) 



we obtain the following mathematical form of the above 
condition for the fractional vapor volume 



1 1 
J /(/') V c (t') + J V c (t',t")df 



df = (27) 



In Eqs. ([24) and ([27]> V c (t', t") is the growth rate dV/dt 
at time t" of a vapor bubble, created with the critical 
size at time t' . Possible effects of "cramming" for many 
neighboring bubbles, not accounted for in Eq. (l27t . are 
not expected to be significant because typically , <k 1 . 
For large values of the Gibbs number 



W c 



» 1 



(28) 
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there is no need to calculate the time integral in Eq. (l27l 
because any integral of the form JT <J>(/')e _G ^' dt' can 
be asymptotically accurately evaluated as 



j <*>(*') 



)e- G{,,) dt' a 0>(f)e 



Vr 



= 0>(f)[-G(f)] V C( " (29) 
by using the first two terms 

G(f') « G(r) + (f - t)G{t) (30) 

of the Taylor series expansion in the exponent — pro- 
vided that the time derivative G = dG/dt < 0. Finally, 
having applied the approximation (l29l to Eq. (1271 1. we 
rewrite our criterion for instantaneous transition from a 
metastable state (p, T) = {p,T\) to the equilibrium state 
(p, T2) in the form 



/3cr 
NVA — 
\mn 



1/2 



d I W c \Y l I W c \ 

&[—) exp r-J 



P/(^2) ~P 

Pi(Ti) - Pv(Ti)' 



(31) 



Note that within the approximation (1291 the bubble 
growth rate V c (t', t") drops out from the answer which 
corroborates the premise that the growth of supercritical 
bubbles can indeed be neglected. 

4.3. Constraints due to applicability limits; the effective 
boiling rate 

Formula (IBTb has been derived under a number of 
physical approximations, some of which break down 
near the spinodal and critical point, as well as on very 
short time scales. As a consequence, to avoid physi- 
cally absurd results, additional constraints may need to 
be imposed on Eq. OTT l if it is to be used over the entire 
metastable region. First of all note that as a (v, p) state 
approaches the spinodal, the difference p e -p grows and 
the critical radius r c falls. Depending on a particular 
model of EOS and surface tension cr (which normally 
decreases with the increasing temperature and decreas- 
ing bubble radius 01), one ma Y formally calculate the 
values of r c less than the intermolecular distance. In 
such a case nothing bad happens with the exponential 
factor in Eq. (1311 , which simply approaches unity, but 
the preexponential factor becomes unphysical: clearly, 
the number of molecules inside a critical bubble, given 
by the product NV C , must at least be larger than one 



more realistic would probably be limiting values NV C ^ 
10-100. 

Less clear is a possible limit of cr — > near the spin- 
odal (whether this should be the general case, still seems 
to be an open question [Q]]). It would hardly be rea- 
sonable to accept that the nucleation rate J, given by 
Eq. (l23l . goes to zero on the spinodal. Then, a satis- 
factory practical choice might be to use in the preexpo- 
nential factors in Eqs. (1231 and OH a fixed (at a given 
pressure p) value of cr = cr m , measured at the melting 
point. The latter implies that formally one can reach 
the spinodal before the MS — > EQ transition occurs; 
in practice, however, because of a very short required 
hydrodynamic time scale, it would mean that the hydro- 
dynamic approach itself is no longer applicable. 

Another possibility is that the effective rate of explo- 
sive boiling 



3cr 
nm 



NV C 
xexp^-j, 



dT\ T I 



(33) 



inferred from Eq. (13H . formally goes to infinity in the 
limit of cr — » for a fixed value of NV C . In this case 
one can use the same model for cr everywhere because 
an infinite value of tJ } 1 along the spinodal and at the 
critical point would be consistent with our principal ap- 
proximation of instantaneous explosive boiling; also, it 
would guarantee that the spinodal itself can never be 
crossed without a transition to the EQ-EOS. 

The effective boiling rate ( 1331 is defined such that for 
the isobaric heating the criterion (|3~H is equivalent to 



d\nT 

dt 



c p T 



(34) 



where dlnT/dt is the full (Lagrangian) time deriva- 
tive in a fixed fluid element. In a general case, when 
the quasi-static approximation is no longer valid, one 
should use the original criterion (|3H . where the time 
derivative in the preexponential factor of Eq. (l3~H is cal- 
culated along the actual evolution path of the considered 
fluid element. 

Two different models of surface tension were used in 
this work to evaluate the effective boiling rate. One is 
the Guggenheim-Katayama law j^L fToh 



o-(D 



T c -T 
T - T 
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(35) 



NV C > 1; 



(32) 



where T c is the critical temperature, T m is the melt- 
ing temperature, and cr m is the surface tension at the 
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Figure 4: Temperature dependence of the effective rate of explosive 
boiling tT 1 , defined in Eq. <33t . for two different formulas for the 
surface tension tr at a fixed pressure p = 500 bars. 



melting point. For SiOa one has T m = 2073 K and 
o~ m = 307 dyne/cm while the FEOS model yields 
T c = 4868 K. 

As an alternative to the Guggenheim-Katayama for- 
mula d35l l. we also used the dependence 



<r{T,p) =<r m 



T c -T\( p,(T)-p v (T) 



T c ~ T m ] \pi(T m ) - p v (T m ) 
P-PI, SP {T) " l/2 



2/3 



Pi(T) - pi, sp (T) 



(36) 



proposed in Ref. 1 12] to make the surface tension vanish 
not only at the critical point but along the entire spin- 
odal; here pi, S p(T) is the density along the liquid branch 
of the spinodal. 

Fig. [4] shows the temperature dependence of the ef- 
fective rate of explosive boiling r^ 1 of SiC>2, calculated 
for the two surface tension models (l35t and (l36t : the 
pressure was fixed at p = 500 bars. Both curves rise 
very steeply over a relatively narrow range of tempera- 
tures to the maximum possible values t7 1 = 10 12 s _1 . 

r b.max 

Intersection with a slowly varying effective heating rate 
T J, 1 ~ l/( c pT) yields the point of explosive boiling 
T = Ti = 3923 K and 3903 K for, respectively, the 
Guggenheim-Katayama law and Eq. ( f36b . This inter- 
section occurs at the value of the Gibbs number G - 15, 
which is large enough to justify the above made approx- 
imations. The difference of 20 K between the two T\ 
values illustrates the accuracy, with which one can cal- 
culate this quantity. 



In the quasi-static isobaric approximation, all the 
fluid elements of a uniformly heated planar foil reach 
the transition temperature T\ and boil up simultaneously 
— which significantly simplifies the analysis of the en- 
suing hydrodynamic picture. In the numerical simula- 
tions (which always exhibit small deviations from the 
quasi-static solution) a simultaneous transition from the 
MS-EOS to the EQ-EOS was simply enforced over the 
entire sample once the boiling criterion was fulfilled in 
its center. The results, discussed below, were obtained 
with the transition temperature T\ = 3923 K, calculated 
with the Guggenheim-Katayama <x law. 

4.4. Subsequent evolution of the two-phase states 

Once a uniform sample undergoes transition to the 
equilibrium state, all its fluid elements start from the 
same point B2 on the p-v diagram of Fig. [3] But then 
the trajectories of different fluid elements may diverge 
because such a transition is typically accompanied by 
a dramatic (by more than a factor of 10) reduction of 
the sound speed c s , and the quasi-static approximation 
may no longer be applicable. As a consequence, the 
temperature and pressure of fluid elements deep inside 
the sample would continue to rise because of persistent 
external heating until the rarefaction wave, propagating 
from the boundary surface, arrives. 

If, however, the condition (flTT i for quasi-static ex- 
pansion is still fulfilled, all the fluid elements would 
quickly (on the time scale t s ) relax to the ambient pres- 
sure p = pq. If this relaxation occurs with a negligi- 
ble contribution from the external heating and without 
shock dissipation, the specific entropy remains constant. 
Thus, the lower limit to the fluid state trajectories be- 
yond point Bo in Fig.[3]corresponds to a quick isentropic 
expansion to point B3, where the pressure pi is equal to 
the initial pressure po and the specific entropy S3 = S2 
is the same as at point B2, followed by a quasi-static 
expansion along the equilibrium isobar p = po. This 
limiting trajectory would also correspond to the fluid el- 
ements at the outer boundary kept at the fixed pressure 
Po- 



5. Hydrodynamics of a planar Si02 foil undergoing 
explosive boiling 

5.1. Numerical setup 

One-dimensional hydrodynamic simulations were 
performed with the CAVEAT code for a planar 
layer of Si02 with the initial thickness of lg = 10 pm. 
The initial thermodynamic state was taken at 7o = 
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Figure 5: Temporal profile of the surface velocity n/(f) of the uni- 
formly heated Si02 foil in two cases: obeying the metastable EOS 
until the explosive boiling occurs at t e t = 52.95 ns (MS case), and 
obeying the equilibrium EOS immediately after crossing the binodal 
(EQ case). 



312 K, po = 500 bar, which corresponds to the ini- 
tial density of po = 2.20 g cirT 3 and the sound speed 
c s = 4.56 x 10 5 cm/s. Starting from t — 0, a constant 
volumetric heating at a rate q = 10 11 J g -1 s was ap- 
plied for the period of < t < 150 ns. Since the initial 
sonic time t s — lo/c s « 2.2 ns is short compared to the 
heating time, the preconditions for quasi-static heating 
are fulfilled. To suppress the unwanted quasi-elastic os- 
cillations, the initial velocity profile was set according 
to the quasi-static solution ©. The foil was simulated 
under a fixed boundary pressure equal to its initial value 
po = 500 bar, which enabled a much clearer visual rep- 
resentation of the boundary fluid-element trajectories in 
Figs. |7]and[8]below. 

The S1O2 EOS was calculated by using the FEOS 



package ( Appendix A i, where the solid part of the ionic 
contribution was omitted for simplicity. Also, the fol- 
lowing input parameters were used: p = 2.2 g cm -3 , 
T = 300 K, and K = 37 GPa (the bulk modulus). The 
parameters of the soft-sphere model, e co i, = 9.74 kJ/g, 
m = 0.8, and n = 2.8, were adjusted such as to make 
the critical point approach the data by Iosilevskiy and 
Gryaznov J3 @]: T c = 6300 K, p c = 3.5 kbar, and 
p c = 0.65 g cirT 3 . 



5.2. Simulation results 

Two cases were simulated: (i) the MS case, where 
SiC>2 obeys the metastable EOS until the temperature 
T\ = 3923 K of explosive boiling is reached at t — t e b = 
52.95 ns, after which a jump to the equilibrium EOS 
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Figure 6: Density profiles in the Si02 foil at four characteristic times 
in the MS case. 



has been made, and (ii) the EQ case, where the con- 
tinuous transition to the EQ-EOS is made immediately 
after crossing the binodal. Fig. [5] shows the evolution 
of the surface velocity ui(t) for both cases. The initial 
oscillations could not be fully suppressed by assigning 
the quasi-static initial velocity profile © because the 
parameter at s had a finite value. 

As might be expected, in the MS case one observes a 
strong jump in the surface velocity ui at the moment of 
explosive boiling due to a sudden increase of pressure 
from po = 0.5 kbar to pz = 1 -6 kbar (the temperature 
jumps by 40 K to T 2 = 3964 K). In the EQ case, on the 
contrary, the surface velocity rises rather gradually on a 
time scale of At 15 ns. Note that, once the artificial 
boundary pressure of po = 500 bar is lowered, the ve- 
locity rise in the EQ case becomes even more gradual 
than shown in Fig. [5] We expect that such a signifi- 
cant difference in temporal behavior of the surface ve- 
locity of volumetrically heated samples should be defi- 
nitely differentiable in corresponding experiments; also, 
it clearly demonstrates the necessity of adequate mod- 
eling of the explosive boiling for interpretation of such 
experiments. 

Fig. [6] presents the evolution of the density profile 
across the foil in the MS case. Prior to the moment 
of explosive boiling t = 52.95 ns the liquid foil is 
bounded by sharp surfaces and has a box-like density 
profile. Having boiled up, it expands as a planar layer 
of gas: had the boundary pressure been set to zero, the 
spatial density profile would have approached the Gaus- 
sian shape, characteristic of the self-similar regime with 
a uniform temperature and linear velocity distributions 
ill- 
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Figure 7: Trajectories of three selected Lagrangian cells of the simu- 
lated SiC>2 foil on the (v, p) phase plane in the MS case. Evolution in 
time is indicated with a series of time markers. 



Figs. [7]and[8]display the trajectories of three selected 
Lagrangian cells on the thermodynamic (v, p) plane as 
calculated, respectively, in the MS and EQ cases. Ini- 
tially, the foil thickness was divided into 100 equal mesh 
intervals, and the simulations were performed in the 
purely Lagrangian mode. Hence, cell 50 corresponds to 
the mid-plane of the foil, cell 25 - to the half-distance 
(along the mass coordinate) between the mid-plane and 
the edge, and cell 3 lies near the foil edge. 

Qualitatively, the v-p trajectories in Fig. Q are simi- 
lar to those depicted schematically in Fig. [3] Up to the 
very moment of explosive boiling, all the fluid elements 
follow roughly the same quasi-isobaric path, but begin 
to diverge sharply after the jump to the equilibrium state 
with p = p2 = 1 .65 kbar. While the liquid-vapor mix- 
ture near the foil edge rapidly expands, with pressure 
relaxing to the boundary value of po = 0.5 kbar, the 
pressure in deep interior of the foil continues to rise up 
to Pmax ~ 4 kbar until the rarefaction wave reaches the 
foil center. 

A salient feature in Figs.|7]and[8]is that, after crossing 
the binodal but before arrival of the rarefaction wave, 
the v-p trajectories of the inner foil elements virtually 
stick to the binodal. Such behavior is caused by a large 
drop in the value of the sound velocity across the bin- 
odal, namely, by a factor ranging from 30 to 3 in the 
relevant pressure interval. Because the sonic time in- 
side the binodal becomes large (tens of nanoseconds), 
one would expect that the pressure in the heated foil 
center should rise practically vertically on the v-p dia- 
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Figure 8: Same as Fig.^]but for the EQ case. 

gram (i.e. at a practically constant density). But then the 
corresponding v-p trajectory would return into the pure 
liquid phase outside the binodal — where, however, the 
sonic time is short (a few nanoseconds) and the foil ma- 
terial would be forced to rapidly reenter the two-phase 
region along a quasi-isobaric path. As a result, for about 
20 ns the binodal becomes sort of an attractor for the v- 
p trajectories of the fluid elements near the foil center. 

Finally note that at t > 100 ns the v-p trajectories 
from the two simulations become quite close to each 
other, i.e. the final state of the evaporated foil does not 
depend strongly on whether one uses the metastable or 
the equilibrium EOS in the two-phase region. It is the 
dynamical behavior inside the two-phase region which 
requires adequate modeling of the boiling (condensa- 
tion) process in the superheated liquid (supercooled va- 
por). 

6. Conclusion 

In this work we investigate a physical model aimed at 
a self-consistent description of the dynamics of matter 
evolving through the two-phase region of the phase dia- 
gram, where the liquid and the vapor states can coexist. 
In particular, we analyze in detail the case of volumet- 
ric heating of a liquid that should eventually result in 
its full vaporization. The proposed model is based on 
the assumption that in the two-phase region the heated 
liquid should be described with the metastable EOS un- 
til the criterion of explosive boiling is fulfilled, at which 
point an instantaneous irreversible transition to the equi- 
librium EOS (obtained by applying the Maxwell rule) 
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can be made. Such an approach appears to be the sim- 
plest possible and still physically self-consistent in the 
framework of pure hydrodynamics, where no additional 
kinetic equations for description of phase transitions are 
to be solved. 

Starting from the theory of homogeneous bubble nu- 
cleation in superheated liquids, we formulate the crite- 
rion for explosive boiling (l3"Tt that has the advantage 
of being local in space and time. We argue that the ir- 
reversible jump to the equilibrium EOS must be made 
under the conditions of constant density and specific in- 
ternal energy, and show that such a jump is accompanied 
by a step-like increase of pressure and temperature. The 
only quantity beside EOS that enters the boiling crite- 
rion OTb is the surface tension of liquid cr, whose de- 
pendence on temperature and density is assumed to be 
known. 

The main properties of the proposed model are illus- 
trated on a simple example of a uniformly heated planar 
layer of SKD2. Hydrodynamic simulations of this layer 
have demonstrated that adequate modeling of the tran- 
sition from the metastable to the equilibrium state can 
be rather important for correct interpretation of experi- 
mental data when, for example, the thermal expansion 
of such foils is measured. 

In the present work the proposed method of hydrody- 
namic treatment of the two-phase region was applied to 
a perfectly uniform sample, where the boiling transition 
occurs simultaneously over the entire volume. However, 
the local transition criterion (|3H can of course be used 
in situations where the distribution of thermodynamic 
parameters is non-uniform (or slightly non-uniform). In 
such cases one may observe a qualitatively different be- 
havior: small initial non-uniformities can be strongly 
enhanced by sudden local jumps in pressure. In reality 
this phenomenon may manifest itself as a spontaneous 
generation of foam-like structures — like it has been 
observed by a quasi-uniform Joule heating of explod- 
ing wires 11511 . But a detailed analysis of such behavior 
remains for future work. 

As a word of caution, it should be noted that, having 
based our main criterion (|3H on the theory of homoge- 
neous nucleation, we have tacitly assumed that possible 
influence of impurities and/or boundary conditions can 
be ignored. Very often this may not be the case, and 
then a more complex analysis, based for example on the 
theory of heterogeneous nucleation may be needed. 
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Appendix A. EOS model FEOS 

The "Frankfurt equation-of-state package (FEOS)" is 
a C++ computer code that can generate EOS tables for 
arbitrary materials. It is an improved version of the code 
MPQeos, which is described in the corresponding MPQ 
Report [16]. The underlying physical model is the "quo- 
tidian equation-of-state model (QEOS)" lfl7ll . A more 
detailed documentation on the FEOS package can be 
found at http://th.physik.uni-frankfurt.de/~faik/ 



The MPQeos model 

The original MPQeos code can only generate EOS 
of chemical elements. Thermodynamic quantities, like 
pressure, specific internal energy, specific entropy, are 
derived from the specific Helmholtz free energy / = 
e - Ts, which is composed of three contributions: the 
electronic part, the ionic part, and the semiempirical 
bonding correction. All three components of / are as- 
sumed to be a function of a single temperature T, al- 
though in principle there is no difficulty in separating 
the ion and the electron temperatures. For any partic- 
ular element, the MPQeos model requires two empiri- 
cal parameters to calibrate the generated EOS, namely, 
the density p„ and the bulk modulus K a = p {dp I dp) 
at a certain reference value of temperature T„ and zero 
pressure p = p a = 0. For substances that are in 
the solid state at normal conditions, usually the values 
T as 300 K are used. 

Electronic contribution. The electronic contribution is 
calculated by using the Thomas-Fermi (TF) model 
iflU, [l9ll . Matter is envisaged as consisting of elemen- 
tary spherical cells, for which the equilibrium electronic 
distribution — and therewith also the thermodynamic 
quantities of the electrons — is calculated by solving the 
TF equation. An important advantage of the TF model 
is the fact that the TF equation and all thermodynamic 
quantities scale with the atomic number Z, and therefore 
must be calculated only once, e.g. for hydrogen. The 
most important drawback of the simple TF model is the 
absence of attractive (bonding) forces between neutral 
atoms. Hence, pressures near the normal conditions are 
overestimated. 
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Bonding correction. MPQeos uses a bonding correc- 
tion from Ref. ll20ll . It adjusts the EOS at a refer- 
ence point (p, T) = (0, T„) to have the correct (or user- 
defined) values of the density p„ and the bulk modulus 
13] K a — p {dp I dp). The reference temperature T is 
usually chosen such that at p = the studied mate- 
rial is in the solid state — i.e. well below the critical 
temperature T c . A relatively small difference between 
the isothermal and isentropic bulk moduli is usually ig- 
nored. The bonding correction has the strongest effect 
on EOS near the solid density. At higher densities the 
TF electronic contribution dominates. 

Ionic contribution. Contribution of ions to thermody- 
namic quantities is fully independent of that from elec- 
trons. It is calculated by using the Cowan model as de- 
scribed in Ref. [Ph . which employs analytical formulas 
to smoothly interpolate between the Debye solid, the 
normal solid and the liquid states. The ionic internal 
energy and pressure are of purely thermal nature. 

Modifications in the FEOS model 

Three important changes were introduced into the 
original MPQeos code, which have lead to the improved 
version FEOS. 

Cold curve improvement. The QEOS model overesti- 
mates the critical pressure p c and the critical temper- 
ature T c . This problem was solved by replacing the 
TF cold curve and the bonding correction for densities 
p < p with a soft-sphere function proposed by Young 
and Corey 121 ] 



e co id (p,T = Q)= Ap" - Bp" 1 + e coh . 



(A.l) 



Here e eo /, is the cohesive energy, or the enthalpy of sub- 
limation per unit mass. The constants A and B are ad- 
justed so as to make the total pressure and the internal 
energy be equal to zero at the reference point (p , T ): 



piPo, T ) = e(p , T a ) = 0. 



(A.2) 



The free parameters m and n are used to improve the 
agreement with the experimentally (or theoretically) 
known critical point. At p — p the sound velocity is 
allowed to be discontinuous. 

Calculation of mixtures of elements. A well known rule 
for the TF mixing of elements fl7ll was added to the 
FEOS. The ionic contribution and the bonding correc- 
tion are handled as a single species with a mean atomic 
number Z and weight A. For the electronic contribution 



of a mixture the partial volumes of all species k are iter- 
atively adjusted in order to equilibrate the TF pressures 
p e ,k and to fulfill an additive volume rule. Finally, the 
appropriate combination of the electronic free energies 
is calculated. 

Calculation of liquid-vapor phase coexistence. In the 
original MPQeos the liquid-vapor phase coexistence 
was calculated for each isotherm below the critical tem- 
perature T c by using Maxwell's geometrical rule (01 . 
Especially for low temperatures T <k T c , this rule is 
computationally very intensive and imprecise. In the 
FEOS this problem was solved by applying a new rou- 
tine, where the saturated vapor pressure is calculated by 
equating the values of the Gibbs' free energy and pres- 
sure on the liquid and vapor branches of the binodal for 
each isotherm below T c . 

Appendix B. Acoustic oscillations of a planar foil 

Consider slow subsonic (u <k c s ) one-dimensional 
motion in a planar foil, which has no heat exchange with 
external sources (q = 0) and whose equilibrium state at 
a fixed boundary pressure po corresponds to a uniform 
density po; the equilibrium foil thickness is Iq. Then 
the density and the pressure perturbations are relatively 
small and can be cast in the form 

p(t,x) ~p +pi(t,x) =p [l +rj(t,x)], 

, , (dp\ . . (B.l) 

p(t,x) * p + ^ — j pi = p Q +p C s T}(t,X), 



where \q\ <k 1. Substituting ( IB. II ) into the hydrody- 
namic Eqs. (0 and ( TT31 and retaining only the first- 
order terms with respect to the small quantities 77 and u, 
we obtain 

dri du 

dt ox 
which yield the classical wave equation 



(B.2) 



8 -±=cf^ (B.3) 



! dx 2 ' dt 2 



dx 2 



dt 2 

for both 77 and u. 

Because the pressure is fixed at the left and right 
boundaries x = ±lo/2, i.e. T](t, +lo/2) = 0, our problem 
is mathematically fully equivalent to the classical prob- 
lem of string vibrations 11221 Ch. V, § 3]. The general 
solution is given by a superposition 

T] = n (t, x) = - V — sin(w„f + </>„) cos — , 

(B.4) 
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Zl ClI„x\ 
u on cos(w„f + 4> n ) sin , (B.5) 

of standing acoustic waves with a discrete set of eigen- 
frequencies 



(i)„ = —?-(2n — 1), n — 1,2, 
h 



(B.6) 



where cf>„ and 

u n tire arbitrary constants. The oscillatory 
solution tB.4) . dB.5l ) is periodic in time, with the period 



t 



2n = 2k^ 2u 



(B.7) 



given by the fundamental frequency <x>\ . 

Now, if we consider the case of a sudden turnon of the 
heating rate q = q m at / = under the condition dTT~b . the 
general first-order solution for the fluid velocity u can 
be obtained as a sum of the quasi-static and oscillation 
terms, 

u(t, x) — ax + u (t, x). (B.8) 

For at s « 1 simple addition of the two linear solutions 
is justified so far as the relative changes of a — Yq m lc 2 
and c s remain small, i.e. at least over a certain time in- 
terval < t <sc a -1 , which may still be large compared 
with the oscillation period t . 

With the initial condition m(0, x) — we must have 
m o (0, +/q/2) = +5«/o and the coefficients u on (at least 



some of them) in Eq. ( IB. 51 ) will be of the order of \u on \ 
al(). By virtue of Eq. (IB .41 , this implies that the spatial 
variation of the density 



alo 

5p - po\T]\ -po — = po (at s ) 



(B.9) 



is one order lower in the small parameter at s — i.e. a 
factor (at s y l » 1 higher — than the spatial density 
variation (fT3T > in the quasi-static solution. Hence, to ob- 
tain a self-consistent first-order solution, we must assign 
a zero initial value r] (Q, x) = for the oscillatory com- 
ponent of the density variation. From this we obtain 
cf> n — for all n = 1,2,.. ., whereas the values of u on are 
uniquely determined by the equation 



U °» Sm [ — J 



(B.10) 



Appendix C. Temperature and entropy jumps by 
explosive boiling of a superheated liq- 
uid 

To determine the direction of temperature change 
by quasi-instantaneous explosive boiling, consider the 
metastable, 6ms (v, T), and the equilibrium, €eq(v,T), 
specific internal energies as functions of the specific 
volume v and the temperature T. Clearly, if the heat 
capacity c v> eq = deEQ/dT is positive and the inequality 



£ms(v, T) > € E q(v, T) 



(CI) 



holds over the entire two-phase region [i.e. within the 
interval v/{T) < v < v,,(T)], then the equation 

cms (v, T) = € E q(v, T 2 ) (C.2) 

implies T2 > T. 

To prove ( 1C.U . consider the function 

<f>(v) = e MS (v,T)-e E Q(v,T) (C.3) 

for some fixed value of T. This function vanishes by 
definition on the binodal, i.e. <p(yi(T)) = <p{v v (T)) = 0. 
If, in addition, we have 



—hr < for all v,(T) < v < v v (T), (C.4) 
dv z 

we necessarily obtain ( IC.U . Since €eq(v, T) is a lin- 
ear function of v for any fixed T < T c , we have 
d 2 €EQ(v, T)/dv 2 = 0, and an inequality 



gWv, T) _ d_ 
dv 2 dv 



< 



(C5) 



will necessarily imply To > T in Eq. ( 1C.2I I. For the van 
der Waals EOS one readily ascertains (1C.51 by calculat- 
ing 

d 2 e(v, T) 6p c v 2 



dv 2 



<0, 



(C6) 



where index c marks the values at the critical point. The 
FEOS used in this work also has this property, as is il- 
lustrated in Fig. IC.9I 

Similarly, if one chooses the specific volume v and 
the specific entropy s as the independent variables and 
ascertains that 



Finally, because at t — we have du„(0,lo/2)/dt = 0, 
the amplitude of the oscillatory term in (IB. 8b is 

u osc = \u (0,l /2)\ = ^alo. (B.ll) 



e M s(v,s)> e E Q(v,s) (C.7) 
over the two-phase region, then the equality 

Cms (v, s) = e EQ (v, s 2 ) (C.8) 
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Figure C.9: Metastable (solid) and equilibrium (dashes) isotherms for 
T = 3700 K on the (v, e) plane for Si02 in the two-phase region as 
calculated with the FEOS model. Validity of inequality jC.lt is clearly 
observed. 



would automatically imply 52 > * because (de/ds) v = 
T > 0. Again, ( IC.71 > can be proven for the van der Waals 
EOS by considering the function 



<p(v) = e M s(v, s) - e E Q(v, s) 



(C9) 



at some fixed value of s. Both <p(v) and its first derivative 



<t> (v) = — 
dv 



Peq(v, s)- p MS (v,s) 



(CIO) 



by definition vanish on the binodal, <p(vi(sj) - 
0'(v/(j)) = 0, 0(v v (j)) = 0'(v v (j)) = 0. Having recalled 
the general thermodynamic identity (01 



df 



(Cll) 



and taking into account that for the van der Waals EOS 
c v ,eq = c v jas (because c v is a function of temperature 
only), one can prove that 



dv 2 



> 



(C.12) 



in the metastable region from the liquid side, i.e. at 
v > vi(s) but before the spinodal. Then, since <f>'(y) 
can cross zero not more than once within the interval 
Vi(s) < v < v v (s), we necessarily have <p(v) > ev- 
erywhere in the two-phase coexistence region. The fact 
that for normal fluids (as contrasted to the so called ret- 
rograde fluids) \>i(s) = 00 (i.e. that normally both the MS 
and the EQ isentropes on the (v, p) plane never actually 
exit the two-phase coexistence region) does not affect 
the above argumentation for the particular case of the 
van der Waals EOS. 
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